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ABSTRACT 

We present the results of near-infrared (NIR) H- and K-band European Southern Observatory SINFONI integral field spectroscopy 
(IFS) of the Seyfert galaxy NGC 1566. We investigate the central kpc of this nearby galaxy, concentrating on excitation conditions, 
morphology, and stellar content. NGC 1566 was selected from our NUGA (-south) sample and is a ringed, spiral galaxy with a 
stellar bar in northsouth direction (PA ~ 5°). The galaxy inhibits a very active Seyfert 1 nucleus but narrow line ratios from optical 
observations in the nuclear region are similar to Seyfert 2 galaxies. The recent strong activity phase, as inferred from strong variablity 
in X-ray to IR wavelengths, makes NGC 1566 an ideal candidate to look for feeding and feedback of a supermassive black hole. We 
present emission and absorption line measurements in the central kpc of NGC 1566. Broad and narrow Bry lines were detected. The 
detection of a broad Bry component is a clear sign of a super-massive black hole in the center. Blackbody emission temperatures 
of ~ 1000 K are indicative of a hot dust component, the torus, in the nuclear region. The molecular hydrogen lines, hydrogen 
recombination lines, and [Fe ii] indicate that the excitation at the center is coming from an AGN. The central region is predominantly 
inhabited by molecular gas, dust, and an old K-M type giant stellar population. The molecular gas and stellar velocity maps both show 
a rotation pattern. The molecular gas velocity field shows a perturbation toward the center that is typical for bars or spiral density 
waves. The molecular gas species of warm F[2(1-0)S(1) and cold ^^CO(3-2) gas trace a nuclear gas disk of about 3" in radius with a 
nuclear spiral reaching toward the nucleus. From the equivalent width of H2(1-0)S(1) a molecular ring with r < 3" can be inferred. 
This spiral seems to be an instrument that allows gas to fall toward the nucleus down to < 50 pc scales. The excitation of molecular 
hydrogen in the nuclear gas disk is not clear but diagnostic diagrams show a distinction between the nuclear region and a < 9 Myr 
old star forming region at the southwestern spiral arm. Possibly shocked gas is detected ^ 2" from the center, which is visible in 
dispersion maps of H2(1-0)S(1) and ^^CO(3-2) and in the 0.87 mm continuum. 

Key words, galaxies: active - galaxies: individual: NGC 1566 - galaxies: ISM - galaxies: kinematics and dynamics - galaxies: 
nuclei - galaxies: star formation - infrared: galaxies 


1. Introduction 

The N uclei of GAlaxies (NUGA) project ( [Garcia-Burillo et aL] 
2003 1 ) started with the IRAM Plateau de Bure Interferometer 
(PdBI) and 30 m single-dish survey of nearby low-luminosity 
active galactic nuclei (LLAGN) in the northern hemisphere. The 
project is ideally suited to map the distribution and dynamics of 
(cool) molecular gas in the inner kpc of LLAGN and to study the 
possible mechanisms for gas fueling at a high angular resolution 
0V5 - 2") and high sensitivity. The ongoing implementation 
of the Atacama Large Millimeter/submillimeter Array (ALMA) 
in the Atacama desert in Chile finally allows the NUGA project 
to expand to the southern hemisphere ( [Combes et al.|2013]|2014| ) 
at an even higher angular resolution (~ 6 - 37 mas, assuming the 
full array is used) and sensitivity. The Spectrograph for INtegral 
Lield Observations in the Near Infrared (SINLONI) adds com¬ 
plementary information to the NUGA goal in the near-infrared 
(NIR, jSmajic et al.||2Q14| ). By maping (hot) gas and the mass 
dominating stellar population we investigate star formation and 
the feeding and feedback of nearby LLAGN in the NIR. 


^ Based on the ESO-VLT proposal ID: 090.B-0657(A) 


1.1. Feeding and feedback in AGN 

Unresolved, powerful, and highly ionizing emission detected in 
the centers of galaxies is thought to stem from accretion events 
onto a supermassive black hole (SMBH). These galaxies are said 
to have an active galactic nucleus (AGN). According to the uni¬ 
fied model of AGN this active nucleus consists of a SMBH sur¬ 
rounded by an accretion disk at scales of up to a few lightdays 
which ionizes the surrounding gas on scales from several light- 
days, in Seyfert galaxies, up to lightyears, in Quasi Stellar Ob¬ 
jects, (e.g., broad line region - BLR) to several hundred parsecs 
(e.g., narrow line region - NLR) and up to even larger scales via 
jets. The torus, a dust and gas mantle, surrounds the AGN on 
parsec to tens of parsecs scales and is thought to be responsi¬ 
ble for Seyfert 1 (torus almost face-on toward the observer) and 
Seyfert 2 (torus almost edge-on) classifications of AGN. The ex¬ 
istence of the torus can be inferred from high column densities 
toward Seyfert 2 AGN and dust blackbody emission in Seyfert 1 
AGN with temperatures up to the sublimation temperature of 
dust (^ 1300 K). 

The host galaxy and its SMBH have been found to exhibit 
tight correlations. These correlations mostly apply to the stel- 
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lar bulge (e.g.,|Magorrian et al.|1998t[Ferrarese & Merritt|2000 


[Marconi & Hunt||20Q3| ). The correlations connect the mass of 

the central SMBH with the mass, luminosity, and kinematics of 
the bulge. Latest studies show that these correlations may de- 
pend on the galaxy classification (e.g., barred galaxies, [Grahamj 
|& Scott|2013| ) and that several galaxies show an over-luminous 
bulge or an under-massive SMBH ( [Busch et al.||2014| ). [Lasker[ 


|et al. ( [2014[ ) show that the total luminosity of the host galaxy 
seems a more robust tracer of the SMBH mass than the bulge 
luminosity. The distribution of the gas in the host is essential 
to understand these correlations, since gas is the progenitor of 
stars which mainly contribute to mass, luminosity and kinemat¬ 
ics of the bulge. Gravitational torques are one of the strongest 
form of force to act on the distribution of the gas on these 
scales. Gravitational mechanisms that exert gravitational torques 
such as galaxy-galaxy interactions (e.g., galaxy merger) or non- 
axisymmetries within the galaxy potential (e.g., spiral density 
waves or stellar bars on large scales) can lead to loss of angu¬ 
lar momentum in the gas. Hydrodynamical mechanisms such as 
shocks and viscosity torques introduced by turbulences in the 
interstellar medium (ISM) can remove angular momentum from 
gas, too. The original NUGA (north) project has already stud¬ 
ied the gaseous distribution in more than ten nearby galaxies 
(~ 4 _ 40 Mpc) with results that show a variety of morphologies 


in nuclear regions, including bars and spirals ([Garcia-Burillo 

et al.[[2005 

Boone et al.[[2007} 1 

Tunt et al.[[2008| [Lindt-Krieg 

et al.pOOS 

Garcfa-Burillo et a . 

2009[), rings ( 

Combes et al. 

2004} [Casasola et al.[[2008 [Com 

bes et al.[[2009 

) and lopsided 

disks (Garcia-Burillo et al. [2003 

Krips et al. 2005} [Casasola[ 


On large scales, hundreds to thousands of parsecs, gravita¬ 
tional torques are the strongest mechanism to successfully trans¬ 
port gas close to the nucleus, whereas visc osity torques can take 
over on smaller scales (< 200 pc, e.g., [Combes et al.[[20()4 


van der Laan et al.|2011[ ). Therefore, large-scale stellar bars are 


citation type (e.g., thermal , non-thermal) of the warm gas and 
the excitation temperature ( [Mouri 11994} [Rodriguez- Ardila et al. 
[2004} [Zuther et al.|2007[ ). The cold gas information can then be 
compared to our results. 

We will use the stellar absorption features (e.g., Sii, CO(6- 
3), Mg I, Nai, CO(2-0)) to get insight into the star formation 
history of the nuclear region (e.g., [Davies et al.|2007| ). Star for¬ 
mation, recent or ongoing, on scales of 0.1 - 1 kpc around the 
nucleus is an important process which is frequently found in all 
types of AGN in contrast to quiescent galaxies (e.g., [Cid Fer- 


[nandes et al.[[2004t [Davies et~ar][2006| [Busch et nr[[2(jl5[ ). The 

debate if outflows from the AGN quench or initiate star forma¬ 
tion is still going on, but it is most probable that outflows can 
show both effects. 


1.2. NGC 1566 

NGC 1566 is a barred, ringed, spiral SAB(s)bc galaxy in the 
Dorado group ( [Mulchaey & Regan|1997} [Reunanen et al.|2002| ) 
harboring a low luminosity AGN (LLAGN) at a redshift of 
z 0.005017 ( [Koribalski et al.|2004] ). The bar is oriented in the 
north-south direc tion (PA~ 5°) and th e ring has a diameter of 1.7 
kpc (see Fig. p]). [Alloin et al.[ ( p~985[ ) detect broad Balmer lines 
with widths of FWHMjj^ = 2400+300 km s“^ They also mea- 


an important feature to transport gas toward the inner Lindblad 
resonance (ILR) (e.g., [Sheth et al.|2005[ ) where the formation of 
nuclear spirals and rings is induced. 

The NUGA sample studies the cold gas distribution. We can 
use the cold gas distribution as a complement to our SINFONI 
NIR observations an d compare it to the sites of hot molecu- 
lar and ionized gas ( [Combes et al.[[MT3} [2014} [Smajic et al.| 


sure an increase in flux of the broad HyS line by a factor of 4 - 5 
within 24 days. They find the narrow line ratio of Htr/HyS to be 
^3.1, which is typical for Seyfert 2 galaxies, however, the broad 
line flux variation and the broad to narrow line ratio Uj3(h)/Uj3(n) 
10 shows that NGC 1566 is in a high activity state typical for 
Seyfert 1.2 galaxies. The light crossing time of the NLR is about 
10^-10^ years, hence, the NLR line ratios suggest that in the past 
few 100 years NGC 1566 has shown more the characteristics of 
a Seyfert 2 galaxy than that of a Seyfert 1. [Kriss et ah] ( [1991[ ) 
measure a redshift of the broad emission lines of 2 00 - 1000 kni 
This might result from gravitational redshift ( Netzer|l977 ) 


2014). This enables us to identify ongoing star formation sites 


and regions ionized by shocks (i.e. super novae (SN) or out¬ 
flows). Furthermore, we are able to compare the distribution of 
cold (e.g., CO, HCN) and hot (e.g., H 2 ) molecular gas that will 
give us a more clear insight on fee ding and feedback o f the AGN 
through its ambient gas reservoir. [Riffel et S] ( [2Q13[ and refer¬ 
ences therein) assume that molecular gas (e.g., H 2 ) traces the 
ambient disk structure and in some cases the streaming motions, 
i.e. the feeding of the SMBH, whereas ionized gas often traces 
outflowing material that is above the galaxy disk plane, i.e. the 
feedback from the SMBH. However, molecular emission lines 
in the NIR (e.g., H2(1-0)S(1) emission) are known shock tracer 
that are found in regions shocke d by jets or outfiows (e.g., [Rifelj 
et al.[[2014[ Davies et al.[[2014[ and references therein). Lately, 


Garcia-Burillo et al.[ ( [2014[ ) fad that cold molecular gas (e.g., 
^^CO(3-2) emission) can be found in outfiows. 

Here, we analyze the interactions of nuclear star forma¬ 
tion sites and the AGN with regard to fueling and feedback, of 
both participants. Using the differently excited H 2 lines (e.g., 
H 2 42.12 yum, 1.957 yum, 2.247 pm) and the hydrogen re¬ 
combination line Bry in K-band and the forbidden transition 
[Fell] 41.644 pm in H-band we are able to constrain the ex- 


but can also be induced by an outflow of optically thick clouds 
at the far side or an infall of clouds at the near side of the BLR 
where the far side is obscured from our view. 

NGC 1566 shows variability from the X-ray to IR wave¬ 
lengths. [Baribaud et al.[ ( |1992[ ) conclude that X-ray flickering 
should occur on time scales of 5x10^ s. They also note that the 
NIR flux variation (heated dust reacting to UV flux changes) 
happens within the central 450 pc from the center at time scales 
similar to the broad Ha variability time scale. The extinction to¬ 
ward the nuclear region of NGC 1566 from NLR and BLR line 
ratios seems to be negligible. [Baribaud et al.[ ( P'9921 ) also calcu¬ 
late a dust mass M^ust ~ 7 x10"'^ Mq and a dust e vaporation ra¬ 
dius of about 47 light-days. [Reunanen et al.[ ( [2002| ) used slit spec¬ 
troscopy on the nucleus of NGC 1566. They detect [Feii] in the 
H-band and the H 2 42.122 pm line in the K-band. They note that 
the lines and the continuum appear weaker along the NLR cone 
than perpendicular to it. They do not detect narrow Bry emission 
only a broad Bry component with an FWHMg^y = 2100 km s“^ 

NGC 1566 is an interesting candidate to investigate the feed¬ 
ing and feedback of the nuclear region and the AGN since the 
galaxy seems to be waking up from a low activity state in the 
past few hundred years (low ionization level in the NLR) and 
shows signs of high activity with strong variability and high ion¬ 
ization on scales smaller than the NLR. [Combes et al.1 ( [2014[ ) find 
from ALMA ^^CO(3-2) and HST data that gravitational torques 
are a very likely cause to drive gas infall in NGC 1566. They 
derive negative torques from 300 pc down to 50 pc. From this 
point dynamical friction can drive the gas to the nucleus. 
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S. Smajic et aL: The nuclear gas disk of NGC 1566 dissected by SINFONI and ALMA 


FOV 10"xl0” 


20” (960 pc) 


Fig. 1 - An HST image of NGC 1566 taken with the WFC3 
UVIS2 F438W filter at a pivot wavelength of 4326 A. North is 
up and east is left. The white rectangle marks our 10" x 10" field 
of view (FOV). 


We will investigate the central 10" x 10" to find if the hot gas 
distribution and kinematics are similar to the cold gas observed 
with ALMA. Are there any signs of a feeding of the SMBH 
(e.g., strong ionization of the nuclear region, streaming gas 
motion)? How strongly is it accreting? Is any strong feedback 
(e.g., jets) visible in the hot gas? Is gas infall accompanied by 
star formation or is star formation rather hindered by the torque 
budget found by |Combes et al.| ( [2014| )? 


This paper is structured as follows: in Sect. 2 we present the 
observation and the data reduction. Section 3 states the results of 
our study and how these were derived. Section 4 discusses the 
results and compares them with literature. In Sect. 5 a summary 
of the results is given and the conclusions we take from our study 
are phrased. 

Throughout the paper we adopt a systemic velocity of 
1504 kms~^ and a distance of 10 Mpc for NGC 1566, following 
[Alloin etaL] ( |1985| ); [Combes et al.| ( |2014| ). 


2. Observation and Data Reduction 


In this paper we present the results of our ESO SINFONI (Eisen- 


|hauer et al.|2003| [Bonnet et al.|2004| ) observation of NGC 1566 
with the Unit Telescope 4 of the Very Large Telescope in Chile. 
The 0'.'25 plate scale with an 8"x8" FOV without adaptive optics 
assistance was used. The average seeing was 0'.'5. To increase 
the FOV to 12" x 12" and minimize the overlap of dead pixels in 
critical areas a dithering sequence of the FOV by +2" was intro¬ 
duced. However, the outer 2" have too low quality data, hence 


for the analysis a FOV of 10"xl0" was used. The dithering was 
done at 9 positions where the central 4" x 4" were observed with 
the full integration time. The gratings used are the H-band grat¬ 
ing at a spectral resolution of R 3000 and the K-band grating at 
a resolution of R 4000. Both bands were observed at a digital 
integration time of 150 seconds and an TST... nodding sequence 
(T: target, S: sky), to increase on-source time. The overall in¬ 
tegration time on the target source in H-band is 2550 seconds 
and in K-band 3000 seconds with an additional 1200 seconds 
in H-band and 1500 seconds in K-band on sky. The G2V star 
HIP 33144 was observed in H-band and in K-band within the re¬ 
spective science target observation. The observing strategy was 
to observe it twice (star in opposite corners of the FOV) with an 
integration time of 2 seconds. The ESO SINFONI pipeline was 
used for data reduction except for atmospheric OH line correc¬ 
tion and final cube creation which were done manually. Detector 
specific corrections were performed manually, as well. For more 
details on the reduction see |Smajic et al.|p014| ). The standard 
star was used to correct for telluric absorption in the atmosphere 
and to perform a flux calibration of the target. A high S/N so¬ 
lar spectrum was used to correct for the black body an d intrinsic 
spectral features of an G2V star ( [Maiolino et al.|199^ . The solar 
spectrum was convolved with a Gaussian to adapt its resolution 
to the resolution of the standard star spectrum. The solar spec¬ 
trum edges were interpolated by a black body with a temperature 
of T = 5800 K. The telluric standard star spectrum was extracted 
by taking the total of all pixels withi n the radius of 3 xFWHMpsF 
of the point spread function (PSF; |Howell|[2000| ), centered on 
the peak of a two-dimensional Gaussian fit. The flux calibration 
of the target source was performed during the telluric correction 
procedure. We referenced the standard star counts at T1.662 yum 
and /I2.159 yum in H- and K-band respectively to the flux given 
by the 2MASS All-sky Point Source Catalogu^ 

To determine the spatial resolution of our observation we 
looked at the radial profiles of the telluric star, the continuum 
emission in H- and K-band, and at the radial profiles of the emis¬ 
sion lines [Fell], narrow Bry, and broad Bry (see Fig.|^. We use 
the broad component of the hydrogen recombination line Bry to 
determine an accurate value of the FWHM of the PSF for the K- 
band because we can measure this value from the science data 
itself. Figure shows that the narrow component of Bry and 
[Fe ii] as well as the telluric stars in H- and K-band show about 
the same spatial extent. From this we infer that the spatial reso¬ 
lution in H- and K-band is similar. We measure a FWHM of the 
broad Bry component of ~ 0'.'59 which corresponds to ~ 29 pc. 
Spatially the PSF shows an elongation in the east-west direction 
slightly rotated by 10° (Fig.[^. 

The correction of the OH lines, the cube alignment, the final 
cube creation, the linemap and spectra extraction was conducted 
using our own IDL routines. 

We use calibrated 350 GHz ALMA data and compare these 
with our NIR results. The final product data is publicly avail¬ 
able from the ALMA archive under project ID 2011.0.00208.S 
(PI: Combes). Details on obser vational setup, calibra tion, imag¬ 
ing and quality can be found in [Combes et al.| ( [2014| ) and in the 
project reports in the ALMA archive. The line cube comprises 50 
channels of 10 km/s width and a beam of 0V64 x 0'.'43 at an PA 
of 123°. The rms achieved is 0.05 mJy beam“^ in the continuum 


^ This publication makes use of data products from the Two Micron All 
Sky Survey, which is a joint project of the University of Massachusetts 
and the Infrared Processing and Analysis Center/California Institute of 
Technology, funded by the National Aeronautics and Space Adminis¬ 
tration and the National Science Foundation. 
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radius [as] 


Fig. 2 - The radial profiles of the telluric star in K- and H-band, the 
K- and H-band galaxy continuum, [Feii], narrow Bry and broad 
Bry. 


and ~ 1.3 mJy beam ^ in the line cubes. Imaging and parts of 
the analysis have been conducted with the CASA software (v3.3 
IMcMuliin et al.|2007D . 

3. Results and discussion 

We resolve the central 485 pc of NGC 1566 at a seeing limited 
spatial resolution of 29 pc using the integral field spectrograph 
SINFONI at the VLT. Narrow and broad components of the hy¬ 
drogen recombination line Bry 42.166 pm and several rovibra¬ 
tional molecular hydrogen lines (e.g. H2(1-0)S(1) 42.12 pm) are 
identified. Additionally, several stellar absorption features in K- 
band (e.g. CO(2-0) 42.29 pm, NaD 42.207 pm, CaT 42.266 pm) 
are detected. The H-band shows a variety of stellar absorption 
features (e.g. Sir 4 1.59 pm, CO(6-3) 41.62 pm) and it harbors 
the [Fell] 41.644 pm line, which is important for NIR line diag¬ 
nostics. The flux and the FWHM of all detected emission lines 
are summarized in tablefor the regions: center r = 1", cen¬ 
ter r = 5", center r = PSF, star formation region (SFr), star 
formation region at PSF sized aperture (SF PSF), and center at 
PSF sized aperture (cPSF), see Fig.j^for more detail. The emis¬ 
sion line FWHMa presented in this paper are corrected for in¬ 
strumental broadening. The cold molecular gas tracer ^^CO(3-2) 
and 0.87 mm emission observed with ALMA are used to com¬ 
pare hot NIR gas and dust with cold gas and dust in the sub-mm. 

3.1. The distribution of gas 

We describe and compare the distribution of the observed gas, 
i.e., ionized [Feii] and Bry gas and molecular H2(1-0)S(1) and 
12 cO( 3-2) gas. 


3.1.1. Ionized gas 

Ionized gas is detected in the [Feii] 41.644 pm, Bry, and 
He I 42.06 pm emission lines in the NIR. 

The forbidden transition [Fe ii] 41.644 pm is slightly blended 
by the CO(7-4) 41.641 pm absorption feature in H-band but is 
very strong. With an FWHM of 350 km s“^ it is the broadest nar¬ 
row emission line at the center (Figs.[^&[3b|). The line map of 
the [Fe ii] emission shows a triangular shape pointing in north. 
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east and southwest directions. The FWHM and equivalent width 
(EW) maps both show elongated features in the east-west direc¬ 
tion. Toward the north the FWHM decreases quickly to about 
100 km s“^ To the south-west it stays for about one arcsecond 
above 300 km s“^ before it drops to ~ 150 km s“^ The eastern 
part extends to about 1'.'5 with an FWHM of about 250 km s“^ 
The [Fe ii] emission is point-like on the nucleus with a stronger 
eastern wing and a weaker south-western wing best seen in the 
EW map (Fig. ^). The EW is 6.8 A on the nucleus and only 
0.6 A in the south-west. A small plateau is detected to the east 
with an EW of ~ 0.9 A. 

The Bry line shows, spatially and spectrally, two compo¬ 
nents. We detect a narrow and a broad Bry emission at the center. 
The broad component shows an FWHM of about 2000 km s“^ 
and is redshifted by 340 km s“^ From the spatial distribution 
of the BLR emission a reliable value for the PSF of ~ 29 pc is 
determined (Sect. |^. The broad component is a clear indicator 
for a SMBH. The narrow component at the center has an FWHM 
of about 230 km s“^ and suggests a region size of 13.5 pc, after 
deconvolution with the PSF width. 

We also detect a spatially resolved strong narrow Bry emis¬ 
sion about 1" southwest from the nucleus (Fig. |^. The off- 
nuclear Bry emission shows a rather elliptical shape in the 
south-east to north-west direction. The FWHM in this region is 
~ 100 km s“^ with a maximum of 150 km s“^ in the south-east 
and a minimum of 70 km s“^ in the north-west. The EW of the 
narrow Bry line at the nucleus is about 1.0 A. At the off-nuclear 
region it is about 2.6 A in the brightest spot and falls to about 
2 A along the ellipse. 

In addition, the Hei 42.06 pm emission line is detected. 
The emission is strongest at the center but is also very promi¬ 
nent in the southwestern region where the strong narrow Bry 
emission is detected. The FWHM peaks on the center with 
250 km s“^ whereas the south-western region shows a broadness 
~ 120 km s“^ similar to that of the Bry. The flux distribution 
and with it the EW of both lines is not as similar. The south¬ 
western emission region shows a peak in Bry in its lower part, 
whereas He i is uniformly distributed over the whole ellipse. The 
EW of He I is up to 0.8 A at the nucleus and up to 1.3 A at the 
southwestern emission region. 


3.1.2. Molecular gas 

The most prominent molecular hydrogen lines that we detect in 
K-band are the H2(1-0)S(1) 42.12 pm and H2(1-0)S(3) 41.96 
pm lines. Other detected molecular hydrogen lines are the H 2 (l- 
0)S(2) 42.03 pm, H2(1-0)S(0) 42.22 pm, H 2 (2-l)S(l) 42.248 
pm, H2(1-0)Q(1) 42.41 pm H2(1-0)Q(3) 42.42 pm lines. 

The molecular emission line H2(1-0)S(1) shows strong emis¬ 
sion on the nucleus and reveals a nuclear spiral structure within 
an r = 3" nuclear disk (Figs. [3^[3T| ). The EW map reveals the 
full extent of the nuclear structure. The arms are clearly visi¬ 
ble, with the eastern arm having a higher EW than the western, 
which is visible in the flux map as well. In EW the center looks 
like a gaseous bar. The spiral does not seem to become narrow 
here, however, this might be a resolution effect. The strong off- 
nuclear Bry emission is situated along the southern edge of the 
western arm, where the EW is lower. Both arms are oriented at 
a PA of ~ 90° and point counter-clockwise. The eastern arm 
and the nuclear region both show the maximum in EW of about 
3.4 A, outside these regions the EW drops to ~ 1.75 A. Both 
arms show concentrated cigar-shaped emission in the parts con¬ 
nected to the nucleus and then turn over by about ~ 90° into 
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Fig. 3 - From left to right: flux [10 W m ^], FWHM (corrected for instrumental broadening) [ km s and EW [A] maps of, from top to 
bottom: [Feii], Hei, H2(1-0)S(1), and narrow Bry emission lines. The regions discussed in this paper are marVprAfti] 
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(a) Flux i2CO(3-2) 


(b) LOSV i2CO(3-2) 



Fig. 4 - Panel |(a)| & |(d)| show the ^^CO(3-2) and con ti nuum at 0.87 mm flux maps in [Jy beam ^ km s and [Jy beam ^], respectively, 
overlayed with H2(1-0)S(1) EW contours. Panels ^ ^ show the ^^CO(3-2) first and second moment maps in [ km s"^] overlayed with 

H2(1-0)S(1) LOSV and FWHM contours, respectively. 


more diffuse emission. The arms do not look like a geometri¬ 
cal spiral close to the center. Inferred from the EW map a bar 
seems to connect the spiral arms toward the nucleus. In the outer 
part the two-arm spiral becomes flocculate and forms a ring-like 
structure at a 2" radius. This ring might corresp ond to an inner 
ILR due to a secondary nuclear bar (Sect. |3.4.^ ) or is rather cre¬ 
ated by gas falling in from the 200 - 300 pc scales toward the 
lower angular momentum transport region at < 150 pc scales 
( [Combes et al.|2014] ). 

The line of sight velocity (LOSV) shows a rotation at a PA 
of ~ 45° that reaches velocities of about +150 km s“^ at a radius 


of about 2'.'5 (for more details see Sect. |3.6.2| ). In addition, there 
is a strong gradient in the central region with the line of nodes 
at a PA ~ 0°. This change in the PA is indicative of a bar or spi¬ 
ral density wave. The FWHM reaches a velocity of 200 km s“^ 
at the center and in regions ~ 2'.'5 to the northeast and to the 
southwest of the nucleus. Along the spiral arms the FWHM is 
-130 km s“^ Along the minor axis of the galactic rotation the 
FWHM does not fall to the width of the spiral arms but stays 
at about 160 km s“^ In the northwest and southeast the FWHM 
drops down to 70 km s“^ (Fig. 3h). 


In the mm-regime we use the ^^CO(3-2) line to compare the 
cold molecular gas distribution to our hot molecular gas distri¬ 
bution derived from ro-vibrational H 2 line emission described 
above. In general, the ^^CO(3-2) emission is very similar to the 
H2(1-0)S(1) emission. In both lines an r = 3" disk with a nu¬ 
clear spiral is detected. The nuclear spiral looks almost identical 
when comparing the ^^CO(3-2) emission and the H2(1-0)S(1) 
EW maps (Fig. and [^. The difference lies in the location 
of the emission maxima. The emission line ^^CO(3-2) peaks at 
connection points of the spiral arms to the center. 

The LOSV-fleld is identical in shape, i.e., maxima and dis¬ 
turbances, and value, i.e., bot h gas es show max/min velocities of 


±150 km s ^ see Figs. 


4b 


& 


13d 


The dispersion of the ^^CO(3- 


2) gas is similar in distribution but differs in value, i.e., the NIR 
H 2 dispersion is lower by ~ 30 km s“^ 

3.2. Gas masses 

From the detected molecular hydrogen lines the warm H 2 gas 
mass can be determined using the luminosity of the H2(1-0)S(1), 
Lh 2 ( 1 - 0 )S( 1)5 and the equation 


H2(1-0)S(3) and all other detected molecular hydrogen emis¬ 
sion lines look very similar to H2(1-0)S(1) in shape and v alue 
(e.g., similar velocities, similar flux distribution, see Figs. |A.2| 
and |A.3| ). Therefore, we use the H2(1-0)S(1) emission line as the 
general description of all molecular hydrogen lines. 


Mh2 = 4.243 X 10“ 


30 / 2.H2(1-0)S(1) 


w 




( 1 ) 


following [Turner & Ostriker|(|1977|>;[S^coville et al.| ( |1982| i; |Wol- 


[niewicz et al.| ( |1998j l; |Riffel et al.| ( |20b8[ >. The warm H 2 gas mass 
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S. Smajic et aL: The nuclear gas disk of NGC 1566 dissected by SINFONI and ALMA 


Table 1 - Emission lines 


Line 

Center r = 1" 

Flux [10-**^ Wm-^] 
Center r = 5" Center r = PSF 

SFr 

SF PSF 

cPSF 

[Fen] 

5.93 ± 0.68 


4.82 + 0.38 

0.21 + 0.09 

0.06 + 0.02 

2.39 + 0.13 

H2(1-0)S(3) 

3.80 + 0.33 

15.15 + 3.18 

2.35 + 0.19 

0.57 + 0.05 

0.14 + 0.01 

0.93 + 0.07 

H2(1-0)S(2) 

1.66 ± 0.18 

5.46 ± 1.35 

1.01+0.10 

0.28 + 0.03 

0.08 + 0.01 

0.39 + 0.04 

He I 

0.51+0.14 


0.42 + 0.09 

0.17 + 0.02 

0.05 + 0.006 

0.23 + 0.04 

H2(2-1)S(3) 



0.12 + 0.07 

0.03 + 0.02 

0.012 + 0.006 

0.05 + 0.03 

H2(1-0)S(1) 

3.41+0.15 

11.21 + 1.2 

2.25 + 0.09 

0.52 + 0.03 

0.13 + 0.008 

0.96 + 0.03 

Bry 

0.71+0.16 


0.60 + 0.10 

0.35 + 0.03 

0.14 + 0.008 

0.31+0.04 

H2(1-0)S(0) 

1.04 + 0.12 

3.68 + 1.06 

0.64 + 0.07 

0.17 + 0.02 

0.05 + 0.005 

0.25 + 0.03 

H2(2-1)S(1) 

0.55 + 0.11 


0.35 + 0.07 

0.11+0.02 

0.03 + 0.005 

0.14 + 0.03 

H2(1-0)Q(1) 

3.68 + 0.35 

33.64 + 8.44 

2.23 + 0.17 

0.74 + 0.11 

0.16 + 0.02 

0.86 + 0.07 

H2(1-0)Q(3) 

3.23 + 0.24 

49.70 + 5.7 

1.80 + 0.12 

0.81+0.11 

0.17 + 0.03 

0.69 + 0.05 

FWHM [ km s-^ 

[Fen] 

321+48 


333 + 34 

130 + 89 

84 + 60 

344 + 24 

H2(1-0)S(3) 

211 + 17 

227 + 43 

211 + 15 

144 + 14 

124 + 12 

207 + 14 

H2(1-0)S(2) 

221 + 28 

214 + 63 

219 + 26 

177 + 25 

183 + 26 

210 + 24 

He I 

150 + 54 


165 + 43 

89 + 19 

82+16 

197 + 37 

H2(2-1)S(3) 



70 + 46 

14+10 

56 + 32 

84 + 44 

H2(1-0)S(1) 

191 + 10 

189 + 25 

199 + 9 

139 + 10 

129 + 10 

204 + 8 

Bry 

165 + 44 


199 + 37 

96+11 

103 + 7 

225 + 31 

H2(1-0)S(0) 

173 + 22 

210 + 67 

178 + 22 

116+16 

112+16 

179 + 25 

H2(2-1)S(1) 

172 


176 

115 

111 

177 

H2(1-0)Q(1) 

166 + 14 

153 + 36 

171 + 12 

115 + 19 

95 + 17 

177 + 12 

H2(1-0)Q(3) 

171 + 14 

150 + 20 

185 + 14 

171+26 

170 + 30 

195 + 14 


Notes. Flux and FWHM table for all regions discussed in this paper (see also Fig.j^for abbreviations). The linewidth of H2(2-1)S(1) was tied to the 
width of emission line H2(1-0)S(0) during the fit, hence no error can be given for the linewidth. The FWHM values are corrected for instrumental 
broadening. 


in a 5" radius aperture, which is all of the warm H 2 in our FOV, 
is derived to 57 M©. To esti mate the cold gas mas s we use the 
conversion factor derived by |Mazzalay et al.| ( |2Q13| ) 

^H2(cold) 




= (0.3- 1.6)x 10^ 


( 2 ) 


(warm) 


We find a cold H 9 gas mass in the central 10"xl0" of (1.7-9.l)x 
10^ M©. Combes et al. ( |2014 ) detect 7 x 10^M© in their r = 18" 
FOV from CO(3-2) observations with ALMA. The values are 
in good agreement since the bulk of the mol ecular mass in NGC 
1566 is located in the inner 6 " (Fig.[^|^and [Combes etal.|2014 
and their Fig. 3). 

The cold gas masses of the central r = 3" disk and the 
cPSF region are estimated from H 2 ( 1 - 0 )S( 1 ) and ^^CO(3-2) 
emission (see Fig. 3g and ^). We measure a ^^CO(3-2) flux 
of ~ 570 Jy km for the central r = 3" gas disk and 
~ 12 Jy km s“^ for the cPSF region. Furthermore, we estimate 
masses of ~ 6.6 x 10^ M© and ~ 1.4x 1 0^ M© respectively, u sing 
the Milky Way conversion values from Bolatto et al.| ( |2013'] ). Us¬ 
ing H 2 ( 1 - 0 )S( 1 ) luminosities we estimate cold gas masses for the 
central r = 3" gas disk and cPSF region of (1.4 - 7.5) x 10^ M© 
and (1.5 - 7.8) x 10^ M© respectively. 


3.3. Emission line regions 

The detection of several narrow emission lines gives us the op¬ 
portunity to analyze the emission at the center of NGC 1566. We 
analyze the ratios of the narrow ionized and molecular emission 
lines with the goal of flnding the nature of their excitation. We 
investigate apertures centered on the nucleus and on the ioniza¬ 
tion region situated ~ 1'.'5 southwest of the nucleus. 


3.3.1. Emission line ratios 


The narrow Bry emission line as well as [Feii] and H 2 ( 1 - 0 )S( 1 ) 
transitions can be used in a diagnostic diagram to disentan¬ 
gle photoionization by young, bright stars and shock ionization 
(e.g., supemovae). Young and bright stars can be found in sys¬ 
tems with recent and strong star formation like starburst galax¬ 
ies. LINER galaxies exhibit high [Feii] and H 2 ( 1 - 0 )S( 1 ) fluxes. 
These species are good shock tracer as they are often found in 
regions of supernovae or outflows/jets. 

We And that the nuclear regions cPSF, r = PSF, and r = 
1" lie on the linear transition relation from SB over AGN to 
LINER (see Fig. |^. The r = PSF region is well situated in the 
AGN regime indicating mixed ionization mechanisms, typical 
for AGN. 

The region SFr lies off the linear correlation seen in Fig. Sit 
is situated in the AGN regime, but under the relation where the 
Seyfert galaxies reside. The ratio log [Fe ii]/PayS puts SFr into the 
Hii galaxy region, however, the ratio logH 2 /Bry shifts it to the 
AGN regime. This behavior can be explained by an H 2 overabun¬ 
dance. As mentioned above, the high H 2 fluxes shift the central 
r - aperture toward the LINER like excitation regime. This 
is the case for SFr as well. The high H 2 flux shifts the SFr region 
from the photo ionization or star formation regime toward the 
mixed ionization or AGN regime. The H 2 over Bry ratio is five 
to ten times higher than in typical H ii galaxies. 

The rich molecular gas disk in the central r = 3" shows 
a variety of H 2 transitions in the NIR (see Tab. [T]). There are 
three main excitation mechanisms for molecular hydrogen in 
the NIR ( |Mouri|19^ [Rodriguez-Ardila et al.|2005| ) which can 
be discriminated with the detected H 2 species: 
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Fig. 5 - Diagnostic diagram of log [Feii]/Pa/5 and log H 2 /Bry for 
the central r = PSF, r = 1" and the star formation region (SFr). 
The conversion factor of 0.744 was used for [Fe ii]/ll .644 pm over 
[Fe ii]/l1.257 pm ( [Nussbaumer & Storey] 1988| ). For the conversion 
of Bry to Pa/5 the case B ratio of 0 .17 was used . Open symbols 
correspond to literatur e values from|Larkin et al.|(|1998| LIN ER), 
|Dale et al.|(|2004| SB),|Rodrfguez-Ardila et al.|{|2004||2005| Syl, 
Sy2). 


i) 


ii) 

hi) 


UV fluorescence (non-thermal) can occur in warm high- 
density gas where highly energetic UV photons from the 


Lyman-Werner band (912-1108 A) are re-emitted by the H 2 
molecules. To distinguish the UV pumping (non-thermal) 
from collisional excitation (thermal) higher level transitions 
need to be detected since the lower levels are populated by 
collisions. 

X-ray heating (thermal) is responsible for H 2 excitation in 
regions with temperatures of < 1000 K. At higher tempera¬ 
tures collisional excitation populates the lower levels. 
Shocks (thermal) can collisionaly populate the electronic 
ground levels of H 2 molecules. The rovibrational transitions 
are populated following a Boltzmann distribution where ki¬ 


netic temperatures can be higher than 2000 K (Draine & 
|McKee|1993T ). 


Both diagnostic diagrams (Fig. show similar behavior 
for the investigated regions. The H2(2-1)S(1)/H2(1-0)S(1) ratio 
ranges from ~ 0.15 to ~ 0.25 for all regions. The central re¬ 
gions, i.e. center r = 1", center r = PSF, and cPSF, show the 
lowest ratio in H2(2-1)S(1)/H2(1-0)S(1) declining with aperture, 
i.e., lower non-thermal UV component. The star forming region, 
SFr, and the PSF sized aperture taken here (SF PSF) exhibit the 
highest value and the ratio is increasing with smaller apertures, 
i.e., higher non-thermal UV component. All regions are situated 
beneath the thermal Boltzmann distribution indicating rather a 
shock than an X-ray origin of the thermally excited molecular 
emission. The line ratio measured for the full FOV (r = 5") is 
not very reliable. This ratio exhibits the highest error bars, prob¬ 
ably due to contamination of the low flux continuum in the outer 
parts of the FOV by OH line correction residuals. 


3.3.2. Level population of the H 2 gas 

The ro-vibrational levels will be populated according to the 
Boltzmann equation when we assume thermal excitation. Then 
the excitation temperature Texc can be derived from the inverse 
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of the gradient of the line fitted to the thermalized levels in the 
graph shown in Fig. [7] These will be on a linear relation if the ex¬ 
citation is thermal. The estimate of the population density can be 
inferred from the observed column density ( [Lester et al.|1988| ): 


^col 


f A An 
Aui hc^ kl' 


( 3 ) 


with fiux / in W m“^, A^i the transition probability 
et al.|199^ , A the rest frame line wavelength, h the 
stant, c the speed of light, and Q the aperture size in radian. In 
thermal equilibrium the ratio of two levels can be written as 


( Wolniewicz 
Planck con- 


N' gj (-AE 


(4) 


with column densities N' and N", statistical weights g'j and g'fl 
Boltzmann constant ks, and the temperature of the thermal equi¬ 
librium T. 

The level population diagram in Fig. [7] compares our mea¬ 
sured level population of the H 2 emission lines to UV-excitation 
models derived by [Davies et al.| ( |2003) ). 

Model 1 is a low density model (nn =10^ cm“^) with cool 
T = 100 K gas and a relatively weak FUV field. In model 2 
density and UV field are increased by one order of magnitude 
and a thermal profile is adopted for the temperature with r^ax = 
1000 K. Model 3 is the same as model 2 but with a maximum 
temperature of r^ax = 2000 K. Model 4 is the same as model 2 
but the FUV field is increased by a factor of 100. And model 5 is 
the high density model with an nn = 10^ cm“^ and temperature 
and FUV field as in model 4. For more details on the models see 
[Davies et al.[ ( [2003[ ). 


3.3.3. Line emission at the nuclear region 

The nuclear line emssion is well situated in the AGN regime in 
the diagnostic diagram in Fig.[^ Interestingly, the r = I" region 
with its slightly bigger aperture moves further toward the LINER 
regime. This trend is caused by aperture effects only. The Bry 
fiux at the center stems from a deconvolved region of 13.5 pc, 
whereas the H 2 emission stems as well from the extended r = 3" 
molecular gas disk. The Bry emission at r = 1" is the fiux in 
the wings of the PSF whereas the H 2 emission is present in the 
central gas disk of up to r < 3" distance from the center. Hence, 
larger apertures will shift the H 2 over Bry ratio toward higher H 2 
fluxes and in this case toward the LINER domain. 

The diagrams in Fig. infer that the central regions move 
with smaller aperture toward the shock model at ~ 2000 K indi¬ 
cating a stronger thermal ionization close to the nucleus. How¬ 
ever, dense gas ionized by UV-fiuorescence can show similar 
emission. The v = 1 transitions are thermalized by collisions 
and with higher density of the gas the v = 2 transitions are 
thermalized as well and hence underpredicted with respect to 
lower density gas excited by UV-fiuorescence (e.g., [Sternberg & 
Dalgarno|1989[[T995} [Sternberg & Neufeld[[ 1999} [Davies et al. 

2003 1 [2005 [ ). Therefore we compare the level population of H 2 
to models from pavies et al.[ ( [2003[ ). The v = 1 transitions in the 
central regions, e.g. cPSF, seem to be thermalized with an ex¬ 
citation temperature of Texc ~ 1800 K (Fig. [^. However, none 
of the apertures taken from the center exhibit values of purely 
thermalized gas. This is shown by the v = 2, J = 3 level which 
lies off the excitation temperature line fitted to the v = 1 lev¬ 
els. This should not be the case for a thermal ionization process. 







































S. Smajic et aL: The nuclear gas disk of NGC 1566 dissected by SINFONI and ALMA 



H2(2-1)S(1)/H2(1-0)S(1) 



H2(2-1)S(1)/H2(1-0)S(1) 


(a) 


(b) 


Fig. 6 - Molecular hydrogen line ratio diagrams. The ratios of H2(2-1)S(1)/H2(1- 0)S( 1) versus H2(1-0)S(2)/H2(1-0)S(0) are shown in (a) 
The ratios of H2(2-1)S(1)/H2(1-0)S(1) versus H2(1-0)S(3)/H2( 1-0)S(1) are shown in|(b)| The curves represent the thermal emission at 1000 - 
3000 K. Vertical stripes represent the region where models by [Black & van Dishoeck| ( |1987| ) predict non-thermal UV excitation. Horizontal 
stripes are thermal UV e xcitation models by |Sternberg & Dalgarno| ( |1989| ). The open, magenta triangle re presents thermal X-ray models by 
[Draine & Woods | ( fl990| ) and t he open, turquoise circle represents a shock model from |Brand et al.|fl989| ). The filled dark bull ets connected 
with a dotted lin e in Fig. are the predicted line ratios from a mixture of thermal and low-density fluorescence models of [Black & van| 
|Dishoeck| ( |1987j ). The first bullet from the left represents 10% non-thermal and 90% thermal UV-fiuorescence, the second 20% non-thermal 
and 80% thermal UV-fiuorescence and so on. 


Hence, other ionization processes, e.g. UV-fluorescence, have to 
be taken into account. The v = 2, / = 5 level seems suppressed 
with regard to the v = 2, / = 3 level and fits the thermal equi¬ 
librium fit. This effect is seen in the X-ray models of [Draine &| 
Woods [ ( [1990[ ) which predict a decrement in the v = 2 / = 5 
level. 

We are not able to disentangle ionization and excitation of 
the gas by the AGN or by stars. 



Energy [10^ K] 

Fig. 7 - Molecular hydrogen level population diagram relative 
to the H2(1-0)S(1) transition. The column density N is given in 
[m"^]. The center r - 1" region marks the energy of the level, 
the other regions are shifted for a better presentatio n in the plot. 
The m odels are for dense UV-excited gas taken from [Davies et al.[ 
P003[ ). The linear fit is to the v = 1 levels of region cPSF which 
yields a kinetic excitation temperature of Tgxc ~ 1800 K. 


3.3.4. Star formation 


The star formation history of the central 10" x 10" of NGC 1566 
does not seem to have been involved in recent star formation due 
to the lack of Hii regions. In fact, Bry emission is only detected 
on the nucleus and in one off nuclear region. The luminosity of 
Bry, Lsry, is proportional to the Lyman continuum flu x and can 
be used as a measure for the star formation rate (SFR) ( [Panuzzo[ 
[et al.|2003}[Valencia-S. et al.|2Q12[ ) in the emitting region 


SFR 


1.585 X 1032W 


Mgyr 


(5) 


The two regions of interest are at the nucleus with an 1" ra¬ 
dius and a luminosity of Lsry = 8.45 x 10^^ W and the lentic¬ 
ular region 1''5 southwest from the center with a luminosity 
of LBry = 4.15 X 10^9 W over an area of ~ 1.33 arcsec^. 
We derive a SFR of ~ 5.3 x 10“^ M© yr“^ at the center and 
SFR ~ 2.6 X 10“^ M© yr“^ at the southwestern region. Addi¬ 
tionally, we can estimate the superno va rate (SNR) in this re gion 
using the [Feii] emission. We follow [Bedregal et ar] ( [2009| ) and 
use two different calibrations 


SNRcal97 = 5.38 ^yr-i (6) 

after [Calzefli| ( [1997[ ) and 

SNRAiH03 = 8.08 ^yr-> (7) 


after [Alonso-Herrero et al.[ ( [2003 [ ). The luminosity of [Feii] 
is measured to be L[Feii] = 7.1 x 10^^ W at the center and 


^[Feii] = 2.51 X 10^^ W at the southwest. The SNRs are ~ 
3.82 X 10“"^ yr“\ and ~ 5.74 x 10“^ yr“^ respectively, at the 
center and ~ 1.35 x 10“^ yr“^ and ~ 2.03 x 10“^ yr“\ respec¬ 
tively, in the southwest. The estimates at the central region are 
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Fig. 8 - Bry contours plotted over the HST image at 4326 A. Note 
the bright emission spots at the SFr emission region at (-1,-1). 


upper limits since star formation and the AGN are responsible 
for the excitation of Bry and [Fe ii] and their respective contribu¬ 
tions can not be distinguished. 

The off-nuclear Bry emission in region SFr is a strong indi¬ 
cator for star formation. The EW of Bry is relatively high here. 
The H2(1-0)S(1) EW in that region goes down with respect to 
the same region in the eastern spiral arm indicating additional 
continuum emission, e.g., young star formation. The diagnostic 
diagram in Eig. [^places this region into the AGN regime, how¬ 
ever, it is off the linear correlation. This is an aperture effect. 
PSE smearing due to the earlier mentioned H ii deficiency at the 
center of NGC 1566 will introduce the shift in the diagnostic 
diagram of this region. The log H 2 /Bry line ratio has values of 
down to -0.2 at the position of the brightest spots in the HST 
images (Eigs.[T^[^. These at least three distinct bright emission 
regions in the HST image are probably the brightest or least at¬ 
tenuated star formation regions. However, the elongated shape 
of the Bry emission there indicates more star formation behind 
the dust and molecular gas of the nuclear spiral. 

In the molecular line ratio diagrams the estimated ratios are 
close to the predicted ratios of [Black & van Dishoeck| ( |1987| ) 
for a mixture of thermal and low-density fluorescence models 
(Eig.[^. The two off nuclear regions, SEr and SE PSE, lie clearly 
off the thermal excitation curve with a tendency towards non- 
thermal UV excitation at smaller apertures. The estimated ratio 
of non-thermal to thermal excitation for the region SEr would 
imply a contribution of ~ 30%. This is a strong hint at young 
star formation taking place at regions SEr and SE PSE. Since 
strong Bry emission is detected in the SE PSE region and the 
[Ee ii] emission is low here, and therefore the SER to SNR ratio 
is high with > 10^, it is probable that the star formation here is 
very young, e.g. < 9 Myr. 

The contribution of non-thermal excitation for the central re¬ 
gions is about 10% - 20% (cPSE, r = PSE and r = 1"). This is 
lower than in regions SE PSE and SEr but still a significant value 
that might hint at star formation at the nucleus. 

The H 2 level population shows that the v = 1 transitions 
scatter around the higher density models (model 2,3,4,5). The 
region SEr tends toward lower density gas, e.g. model 1, and 
indicate s the characteristic ortho-to- para shift of fluorescent ex¬ 
citation ( [Sternberg & Neufeld|1999| ). The v = 2 transitions show 
a similar effect but the differences in the models are here more 
evident than for the v = 1 transitions. 
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3.4. Continuum 

We analyze the emission of the 0.87 mm continuum observed 
with ALMA and the NIR continuum observed with SINEONI. 


3.4.1. The millimeter continuum 


The mm-continuum at 0.87 mm (Eig.|^ peaks in the same re¬ 
gion as the NIR continuum but the distribution is different com¬ 
pared to the NIR (see Sect. [3.4.2| ). The 0.87 mm emission is dis¬ 
tributed similar to the molecular lines. Apart from the peak in the 
center, 0.87 mm emission is detected in the south-western spiral 
arm at the position of narrow Bry emission and in the north¬ 
east coinciding with the regions that show an increased width 
in the molecular lines (Eig. |^. [Combes et al.[ ( [2014[ ) And that 
the 0.87 mm continuum is dominated by dust emission in the 
full 18" X 18" EOV. To decide on the dominating mechanism 
in local emission regions is not possible due to the lack of high 
resolution radio data at other frequencies. 


3.4.2. The NIR continuum 

In the NIR the H- and the K-band emission was observed. The 
continuum flux density is stronger in H- than in K-band. Towards 
the center the continuum becomes redder, as is expected in a 
Seyfert 1 galaxy, hence the slope becomes flatter but it is not 
inverted. The H-K map (Eig. 0 shows a clear reddening toward 
the center with an H-K value of > 0.8 mag at the very center. 
This implies that we see warm to hot dust emission in the galaxy 
center ( [Eischer et al.[[2006l [Busch et"ar][2014[ ). Since we see a 
broad Bry component we assume that we as well see the inner 
edge of the dust torus surrounding the AGN. 



2 1 0-1-2 

A RA [as] 


Fig. 9 - The H-K color diagram in magnitudes of the central 4" X 
4" of NGC 1566. The red ellipse at (1'.'5, -1'.'5) shows the beam 


size. For details see Sect. 3.4 


To analyze the NIR continuum further a decomposition was 
performed. The continuum components: stellar template, hot 
dust blackbody, power-law and an overlaying extinction com¬ 
ponent were used to determine the continuum composition. The 
stellar co mponent was f itted u sing resolution adapted template 
stars from [Winge et al.[ ( [2009[ ). Since these stars are only avail¬ 
able in K-band from 2.2 pm to 2.4 pm the decomposition had to 
be performed in this wavelength range. For more details on the 
decomposition see [Smajic et al.[ ( [T012[ ) and [Smajic et al.[ ( [2Q14 l. 

The decomposition provides two interesting results, see 
Fig. [iq[a)|(c)|(f)[ The blackbody, responsible for fitting hot dust, 
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(a) Center r - I" (b) Center r = 5" 
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(e) SF PSF 
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(c) Center r = PSF 
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Fig. 10 - The continuum decomposition for six different regions. The axes are flux density [10 W m ^ /rm versus wavelength [jim]. 
The components are: Red curve’, hot dust component; Yellow curve: power-law component; Full blue curve: total stellar component; Lined 
blue curve: MOIII star; Line-dot blue curve: K3III star; Triple-dot-lined blue curve: G5III star. The combined components curve is in green 
and the original input spectrum is in black. The dotted black curve around zero is the residual of spectrum minus combined curve. Note that 
the emission lines at /12.223 pm and X2.2A1 pm were masked. 


returns a hot ~ 1000 K component with a significant contribu¬ 
tion at the very center. The stellar component produces best re¬ 
sults when HD2490, an MOIII star, representing the old stellar 
population, and HD1737, a G5III star are used. The contribution 
of G5III-like stellar emission becomes apparent at the center. 

In the NIR, the continuum emission of NGC 1566 in the cen¬ 
tral 10" X 10" is mainly produced by stars. The continuum de¬ 
composition, however, reveals a hot d ust co mponent at t he nu - 
cleus. The fits at the center, r = 1" (Fig. |10a| ), r = PSF (Fig. |10c| ), 
and center PSF (cPSF) (Fig. |10f| ), need a hot dust component. 
The hot dust exhibits only 4% of the flux density in the r = 1" re¬ 
gion, whereas the r = PSF region exhibits 7% of the flux density 
and the cPSF region exhibits even 16% flux density as hot dust 
emission at a temperature of ~ 1000 K. The hot dust emission 
is not strong enough to create a steep red continuum in K-band. 
The difference in percentage here is caused by the aperture size 
of the central regions. All regions exhibit a hot dust flux density 
of ~ 500 X 10“^^ W m“^ pm~^ (Fig. 10). This hot dust emission 
at the nucleus is visible in the H-K map as well (Fig. [^. The 
nucleus of NGC 1566 shows characteristics of Seyfert 1 nuclei 
on smaller scales , i.e. broad hydrogen lines, and hot dust con¬ 
tinuum emission. The effect of aperture on the features is essen¬ 
tial. Apertures of more than 100 pc or 2" would probably not be 
able to detect and measure the hot dust component. NGC 1566 
would then show the characteristics of a quiescent galaxy with¬ 
out any narrow H ii emission at the nuclear region, the broad Bry 
emission might be detected further out, depending on its strength 
(e.g., [R^nanen et al.|2002| ). 


The stellar continuum in our FOV is best fitted by an MOIII 
star, a giant. The central r ~ 1" needs a G5III giant as well 
for the fit. The G-type star contributes one third of the flux in 
the central r - V' region. Outside this region the G-type star is 
not needed. The differences between the two stars are the tem¬ 
perature, i.e. the spectral slope, and the depth of the absorption 
features, e.g. CO(2-0). The G5III stellar contribution indicates a 
need for bluer continuum emission or for shallower features in 
the continuum emission at the center of NGC 1566. The need for 
bluer continuum emission might stem from uncertainties in the 
hot dust contribution, however, the hot dust flux stays constant 
over several apertures. The main argument for a G-type star here 
is the EW of the CO absorption feature. The EW in CO(2-0) of 
the MOIII star is too high, i.e. the CO absorption feature is too 
deep, although a hot dust contribution is taken into account. A 
higher hot dust contribution can be excluded since the spectral 
shape of the residual galaxy spectrum, after subtraction of the 
hot dust black body emission, becomes too blue for an K- or M- 
type star to fit. Star formation at the center of Seyfert galaxies 
and in the torus itself is not unusual (e.g. [Davies et al.|2007| ) and 
can account for additional featureless continuum. 

The spatial shape of the continuum emission is roundish with 
a slight elongation at a PA of 34°. This is the same PA as the 
line-of-nodes of the stellar velocity field. Hence, the stars in the 
center do not follow the main stellar bar distribution, which is 
at a PA of ~ 0°. The 34° that we measure are either created by 
a small nuclear bar within the r = 9".l pseudo-ring (Comeron 
jet al.|2010| ) or the angle stems from the bulge population and is 
affected by projection effects. However, to disentangle this a so- 
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(a) K-band continuum 


(b) Stellar continuum 


(c) Stellar dispersion 


Fig. 11 - From left to right: K-band continuum [10 W m ^ pm ^], the fitted stellar continuum [arbitrary units] and the SVD [ km s 


obtained from the continuum decomposition. For more details, see Sects. 3.4 


phisticated decomposition of the galaxy structure using high res¬ 
olution imaging in the NIR is required (e.g. |Busch et al.|2()T^ . 

The stellar velocity dispersion, LOSV, and intensity distribu¬ 
tion are shown in Figs.[^[^ and[^ The stellar LOSV shows 
a rather regular rotational field with a PA of 214° at velocities 
of +60 km s“^ with a boxy redshifted side. The stellar velocity 
dispersion shows velocities of 90 to 125 km s“^ with the very 
center being at 100 km s“^ and does not show signs of kinemat¬ 
ically decoupled regions. 



Fig. 12 - Stellar velocity dispersion histogram. The X-axis shows 
the fitted dispersion in 5 km s~^ bins. The Y-axis shows the total 
number of spatial pixel that correspond to the dispersion bin. The 
red curve is a Gaussian fit to the distribution. 


3.5. The active nucleus 

The AGN of NGC 1566 is a Seyfert 1 nucleus that exhibits broad 
hydrogen emission lines and is showing variability across the 
whole wavelength range. 


3.5.1. Mass of the SMBH 


The broad Bry line flux and FWHM are used to estimate the 
mass of the central SMBH of NGC 1566. A conversion factor of 
Paof/Bry ~ 12 is used to translate the Bry to Pacr flux and be able 
to use the equation 


M. = 10 


J.29+0. 


•1 / 


\ 10"^^ erg s“^ 


0.43+0.03 


FWHMp^, 
lO^km s“^ 


1.92+0.18 

Mo 
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( 8 ) 


derived by |Kim et aL] ( |2010| ). We estimate a black hole mass 
of M, = (3.0 ± 0.9) X 10^ Mq using a luminosity of Lgry = 
5 X 10^^ erg s“^ and a FWHMp^^y = 2000 km s“^ This value 
is less than a factor two smaller than the 5 x 10^ Mq estimated 
from optical broad line measurements ( [Kriss et ar]|1991| ). This 
difference might be introduced by variations in the activity of 
the AGN. A flux increase in the broad HyS line was measured by 
a factor of four to five within only 24 days ( [Alloin et al.|1985| ). 
Assuming a low activity state during our observation, an increase 
of a factor four in luminosity will increase the BH mass derived 
from the Bry broad line to (5.5 ± 1.7) x 10^ M©, a factor of 
two. Hence, the derived value of (3.0 ± 0.9) x 10^ M© is in good 
agreement with the literature value. 

To use the M-cr^ relation and estimate the mass of the 
SMBH the dispersion of the bulge is needed. To determine the 
bulge dispersion a Gaussian fit was performed on the distribution 
of the stellar velocity dispersion values (Fig.p^. The flt yields a 


dispersion of 105 ± 10 km s ^ Following Gultekin et al. ( 2009| 


M. = 10^- X 


200 km s ^ 


4.24+0.41 


M©, 


(9) 


the mass of the SMBH is then estimated to M, = 8.6 + 4.4 x 
10^ M©. This value is similar to the 8.3 x 10^ M© from | Woo & 


|Urry| ( |2002| ) but a factor two to three higher than the BH mass 
estimate from the broad emission lines. 

Graham & Scott] ( |2013| ) have investigated the M-cr* rela¬ 


tion for barred galaxies. They And that barred galaxies follow a 
slightly different M-cr,. relation than non-barred galaxies. Their 
best flt 


M. = X 


cr He 


200 km s ^ 


5.29+1.47 


M©, 


( 10 ) 


returns an upper limit for the BH mass of M, = 6 x 10^ M©. 
When we use their flt with the least root-mean-square scatter 


M. = 10 


J.78+0.1 


200 km s 


-1 


4.14+0.55 


M©, 


( 11 ) 


the result is a BH mass of M, = (4.2 ± 2.4) x 10^ M©. The latter 
two results agree within the errors with mass estimates using the 
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I 



(a) Stellar LOSV 


2 0 

A RA [as] 

(b) Model LOSV 



(c) Residual LOSV 



A RA [as] 

(d) H2(1-0)S(1) LOSV 


(e) Model LOSV 


(f) Residual LOSV 


Fig. 13 - Top panel: Stellar LOSV, fitted LOSV model and the subtracted residual in units of [ km s"^]. Bottom: H2(1-0)S(1) LOSV, fitted 
LOSV model and the subtracted residual overlayed with contours of the EW map of H2(1-0)S(1) in units of [ km s"^]. The EW map was 
overlayed to display the nuclear spiral onto the residual velocity field of the gas. The straight lines denote the orientation of the major and 
minor rotation axis as determined by the LOSV model. Eor more details, see Sect. 13. 6 


broad emission lin es, i.e., broad Bry (see Sect. |3.1.f] ), broad Ha 
( [Kriss et al.|1991| ). 

The scatter of the individually estimated BH masses is higher 
than their uncertainties. Therefore, we derive a mean BH mass 
from the four relations used in this paper to estimate the BH mass 
of NGC 1566. The mean BH mass is (5.3 + 2.9) x lO^Mo. All 
individual BH mass values are within the standard deviation of 
the mean BH mass. 


3.5.2. Activity of the SMBH 


NGC 1566 is actively accreting mass as can be inferred from 
the observed variabilities. The 2-10 keV X-ray luminosity of 
the AGN was measured by ILevenson et al.l (|2009|) to Lx,d 20 = 


1041.5 g -1 Pqj. ^ (Lstance of 21.2 Mpc, hence we estimate an 

Lx = 7. Ox 10"^^ erg s“^ for 10 Mpc. From the derived mean BH 
mass an Eddington luminosity of LEdd = (6.6+3.7)x 10"^^ erg s“^ 
can be determined. Using the rela tion Lhni ^ 16 x Lx = (1.1 ± 


0.4) X 10^2 erg s'^ for LLAGN (Ho 2008 2009) the Edding¬ 


ton ratio of the active nucleus in NGC 1566 can be estimated 
to AEdd - (2 + 1) X 10“^. Thi s is a typical value for S eyfert 1 
LLAGN ( Ho||2008| ). Although |Kawarnuro et al. ( 2013| ) use dif¬ 
ferent parameters, e.g. BH mass, distance, Lx conversion factor 
etc., they find a similar value of 3.2 x 10“^. 


The mass accretion onto the SMBH can then be estimated 
with 


&M 7y|3ol 


dt 


TjC 


2 ’ 


( 12 ) 


where rj is an efficiency factor which is usually of the order 
of 0.1. We calculate a mass accretion rate of (2.0 ± 0.7) x 
10“"^ M© yr“^ using L^oi from above. The cold H 2 gas mass in 
the 10" X 10" EOV was measured to be ~ 5.4 x 10^ M©. This 
mass is enclosed in the central r = 3" gas disk (see Eig. [3^ 
which then has a column density of ~ 1.9 x 10^ M© arcsec^ 
or ~ 7.7 X 10^ M© pc“2. We derive a gas mass density of 
3.6 X 10^ M© pc“^ from the H 2 cPSE flux measurement by as¬ 
suming that the emission stems from a spherical region 29 pc in 
diameter. Hence, the black hole of NGC 1566 has enough mass 
enclosed in its central few parcsecs to accrete at the current rate 
for a few Myrs. 

Due to the narrow line ratio of HajH/S ~ 3.1 ±0.3, which 
are typical for Seyfert 2 galaxies, [Alloin et ar] ( |1985| ) conclude 
that NGC 1566 might be a waking up Seyfert 1, which was in 
a low ionization state in the past few hundred years as inferred 
from emission line ratios in the NLR, beginning the harder ion¬ 
ization of the nuclear region. The AGN would have enough fuel 
to do so and start ionizing hydrogen outside the ~ 13.5 pc at the 
center (see Sect. 3.1.1|). H owever, on the one hand, the line ratio 
from [Alloin et aL ( 1985| ) might have been created by an aper- 


Article number, page 13 of 19 







































































A&A proofs: manuscript no. master 


ture effect and the H ii deficiency in the central region of NGC 
1566 ( [Comte & Duquenn'o^|1982| ). On the other hand, the Hii 
deficiency might result from the low activity state in the past 
of NGC 1566 since there is hydrogen in the nuclear region but 
in molecular form mainly. The activity of NGC 1566 over the 
last ~ 20 yr is a variable one (e.g., [Alloin et al.||1985t [Krissj 


jet al.|199r| [Baribaud et al.|1992l [Levenson et al. 12009 j ) without 
a specific trend of in- or decrease. Hence, we can confirm that 
NGC 1566 hosts a Seyfert 1 nucleus at its center and it has a gas 
reservoir in the central few pc to keep the central engine fueled 
for several Myrs. 

3.6. Kinematics 

The kinematics in the center of NGC 1566 show rotations at a 
PA of ~ 214° for the stellar kinematics (Fig. jlSaj ) and at a PA 
of ~ 222.5° for the molecular gas (e.g., H2(1-0)S(1), Fig. jlSd ). 
The observed velocity field of the molecular gas shows strong 
deviations from rotation at the center. 

The orientation of the galaxy toward the observer can be in¬ 
ferred from considerations of the winding sense of the spiral 
arms (Fig. [T]). These kind of spiral arms are always seen to be 
trailing. Hence, the near side has to be in the northwest when 
the arms are trailing because there is blueshifted motion in the 
northeast and the arms are oriented as seen in Fig.[^ 


3.6.1. Stellar kinematics 

The line of nodes of the stellar rotation is aligned with the stellar 
continuum major axis at a PA of 34°. We derive a PA from the 


H2(1-0)S(1) line of 42?5 in good agreement with jAguero et al. 
( |2004| ). They find a PA of 44° + 8° from optical emission lines. 
The stellar LOSV map shows a smooth rotation field. We fit a 
model to the observed velocity using the Plummer potential to 
represent the bulge gravitational potential ( [Barbosa et al.|2006| ). 
The model subtracted map shows low residuals (Fig. [13c[ ). The 
stellar velocity dispersion shows a slight drop at the very cen¬ 
ter ( Sect. ja. This is often an indication for a ste llar disk 
(e.g. [Emsellem et al.|2Q01t[Falc6n-Barroso et al.|2QQ6[ ) at scales 
< 100 pc. However, there are no features in the stellar residual 
LOSV map (Fig. [13c[ ) that would support a stellar disk at these 
scales. 

The difference in stellar and gaseous PA might be explained 
by streaming motions of the gas related to a strong spiral wave. 
But the gaseous PA over the central 8" x 8" is in agreement 
with measurements at larger scales, hence the effect has to stem 
from larger scale spiral density waves rather than the nuclear 
spiral discussed here. To substantiate any misalignment between 
gaseous and stellar disk, however, a much better estimate of the 
PAs is needed. Hence, we find that gaseous and stellar disk are 
oriented very similar but might be misaligned. 

3.6.2. Gas kinematics 


On scales of 200 - 300 pc [Combes et al.[ ( [2014[ ) showed that 
the gravitational torques are able to transport almost half of the 
angular momentum of the gas in one rotation period. Here gas 
is able to move in spiral arms that connect the inner Lindblad 
resonance (ILR) of the nuclear bar with the central < 200 pc. 
On smaller scales the angular momentum transport is smaller 
per rotation period, but the period becomes smaller as well. The 
inner 200 pc show the molecular nuclear gas disk which exhibits 
a clear two-arm spiral structure (Figs.|^[^ and[4aj). 


The nuclear spiral disturbs the velocity field of the molec¬ 
ular gas and creates a strong S-shaped feature at the center. To 
remove the rotational part from the spiral disturbance we again 
fit a model to the observed velocity using the Plummer poten¬ 
tial to represent the bulge gravitational potential. The residual of 
the fit (Fig. [13f[ ) highlights the difference in the observed LOSV 
(Lig. [13d[ ) and the model (Lig. [13e[ ). The center of the resid¬ 
ual map shows blueshifted residuals to the west and redshifted 
residuals to the east of the nuceleus. The nucleus is situated at 
0 km s~^ Contrary to the case of NGC 1068, discussed by 
[Garcia-Burillo et aT| ( [2014[ ), where non-virial outflow motions 
are required to explain the residual velocity field derived for the 
nuclear 200 pc gas disk of that Seyfert 2 galaxy, in NGC 1566 
residuals shown in Lig. [13f[ can be accounted for by streaming 
motions produced by the spiral/bar structure. The S-shape shows 
non-circular motion in one surface as it is the case for bars or 
warps. Non-circular orbits, e.g. a closed elliptical orbit with axes 
not parallel to one of the symmetry axis (minor or major), can 
produce residual velocities as observed in Lig. [13f[ 

There are no signs for an outflow from the center. [Schmitt[ 
[& Kinne}^ ( [1996[ ) detect a faint one-sided cone in [O iii] pointing 
toward the southeast. The cone is smaller than 0'.'5 and can prob¬ 
ably be associated with the NLR of NGC 1566. The increase of 
the LWHM of H 2 along the minor axis is in the same direction as 
the [O III] cone but is probably caused by beam smearing along 
the 0 km s“^ velocity gradient. At a size of r ~ 0'.'5 the [Oiii] 
cone occupies spatially the nuclear Bry and [Le ii] emission re¬ 
gion, i.e. these lines originate from the same region. However, 
we do not detect one-sided emission on the nucleus from any 
emission line (Lig.|^. The [Le ii] emission has a triangular shape 
but is centered on the nucleus. There is extended flux toward the 
southeast and also a broadening of the [Le ii] line. These are only 
hints at a possible weak outflow from the nuclear region of NGC 
1566 and can as well be associated with the NLR. The LOSV 


residuals up to 1" from the center in Lig. 13f fit an outflow along 
the minor axis, but these residuals are easily explained due to 
deviations introduced by the density waves of the nuclear spiral. 
Higher angular resolution observations are needed to be able to 
compare the [Oiii] emission with shock tracers in the NIR (e.g., 
H2(1-0)S(1) and [Leii]) and look for signs of outflowing gas in 
the LOSV and LWHM. 

At 1'.'5 east and 1" north the dispersion increases over the 
nuclear disk average in both investigated molecular species, i.e. 
H2(1-0)S(1) and ^^CO(3-2), see Lig.|^ No increase in line flux 
can be measured here, but the 0.87 mm continuum shows some 
substantial emission. [Combes et al.[ ( [2QllT ) find that in their full 
18" X 18" LOV the 0.87 mm continuum is dominated by dust. 
Hence, it is most probable that heated dust, rather than syn¬ 
chrotron emission from supernovae or free-free emission, is cre¬ 
ating the 0.87 mm continuum emission. In combination with the 
higher dispersion in these regions it must be turbulences in the 
gas that warm up the dust. What is creating this turbulence is not 
clear. As mentioned above an outflow from the nucleus is rather 
improbable as the only signs for an outflow are hinting toward 
the southeast. Other possibilities are interactions of nuclear disk 
and spiral arms, e.g. the nuclear spiral. However, these are only 
speculations since our data is not giving any substantial hints 
toward the origin of this turbulence. 


4. Conclusion and summary 

We have analyzed NIR ILS data of the central 10" x 10" of 
the Seyfert 1 galaxy NGC 1566. The Hii deficiency reported by 
Comte & Duquenn^ ( [1982[ ) is confirmed. We make a first de- 
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tection of narrow Bry emission at the center and at one region 
offset to the southwest, in contrast to [Reunanen et al.| ( |2002| ). 
They probably were not able to detect the narrow Bry emission 
due to the width of their aperture since the slit positions do coin¬ 
cide with the positions of narrow Bry emission. 

From the detection of a broad Bry component we estimate 
a BH mass of (3.0 + 0.9) x 10^ Mq similar to othe r BH mass 
measurements from broad lines ( |Kriss et al.||199l| ). From our 
continuum decomposition we derive the mass dominating stellar 
distribution and its velocities. Using the velocity dispersion of 
the stars in the bulge we estimate the BH mass independently 
from broad emission lines. We find a BH mass of (4.2 + 2.4) x 
10 ^ Mo, in excellent agreement with the BH estimate from broad 
emission lines. We find a BH mass of (5.3 + 2.9) x 10^ Mq as 
a mean value from all methods and relations used in this paper. 
Furthermore, the Seyfert 1 classification is verified due to hot 
dust blackbody emission of ~ 1000 K at the nuclear position, 
which indicates that the hot inner edge of the torus structure is 
visible. 

The accretion onto the SMBH is typical for LLAGN with Ed¬ 
dington ratios of AEdd = (2 + 1) x 10“^. The estimated accretion 
rate is ~ 10“"^ and therefore 10^ to 10"^ times higher than for 
SgrA* ( [Bower et al.|2003[|Nayakshin|2005| ), which has a simi¬ 
lar SMBH mass. 


Molecular hydrogen is very strong in the observed FOV par¬ 
ticularly in an r = 3" disk and in a spiral structure within this 
disk. Disk and spiral are both detected as well in the warm H 2 
emission lines with SINFONI as in cold H 2 gas from the ^^CO(3- 
2) emission with ALMA. The shape of the spiral looks similar 
when comparing the EW of H2(1-0)S(1) with the flux map of 
^^CO(3-2). However, H2(1-0)S(1) shows the strongest emission 
on the nucleus whereas ^^CO(3-2) shows several strong emission 
spots. The southwestern emission spot of ^^CO(3-2) coincides 
with the Bry emission seen in the SF region. From our H 2 emis¬ 
sion we determine a cold H 2 gas mass of (1.7 - 9.1) x 10^ Mq 
which is in agreement with results based on the ALMA observa¬ 
tion by [Combes et al.| ( [2014| ). The molecular gas seems to form 
a ring-like structure at r ~ 2" best seen in the EW map in Fig.[3T| 
and the ^^CO(3-2) emission line map (see Fig. 4a). 

Along the southern spiral arm a star forming region is de¬ 
tected at a distance of ~ 1'.'5 from the center. The measured SFR 
over a surface area of ~ 1.33 arcsec^ is 2.6x 10“^ Mq yr“^ Com¬ 
paring this to the global Schmidt law ( [Kennicutt|1998[ ) we find 
that our data point for SFr lies very close to the relation. How¬ 
ever, the value for the nuclear 3" x 3" disk is situated below the 
relation with an SFR of 8.0 x 10“^ M© yr“^ over a surface area 
of TT X 9 arcsec^. The nuclear region of NGC 1566 has a large 
molecular gas reservoir which it is not using efficiently for star 
formation. 


The excitation mechanism of the molecular gas is partly due 
to thermal processes (e.g., shocks) and partly due to non-thermal 
excitation (e.g., UV fluorescence). The diagnostic diagrams in 
Figs. [^&|^ show a clear distinction between the nuclear region 
and the star forming region. The apertures taken from the nu¬ 
clear region are situated in the AGN regime in the diagnostic 
diagram in Fig.|^ The SF region is situated right of the starburst 
region due to the above mentioned H 2 overabundance. The H 2 
line ratios (Fig.[^ imply young star formation at regions SFr and 
SF PSF due to a non-thermal excitation fraction of up to 30%. 
From the strong Bry but very low [Fe ii] emission a young star- 
burst of < 9 Myr is suggested. Therefore we find that gas infall 
is accompanied by young star formation that can be associated 
with the nuclear spiral arms. 


The central regions are closer to the shock models but show 
a non-thermal component as well (Fig.[^. The level population 
diagram (Fig. [ 7 ]) shows that the v = 2 levels are not thermalized 
with the V = 1 levels indicating that the region is excited by a 
combination of thermal and non-thermal processes in dense gas, 
i.e. shocked gas and gas ionized by UV-fiuorescence from young 
newly formed bright stars. Both methods hint on star formation 
at the nucleus which we are not able to distinguish further from 
ionizing emission of the AGN. 

A region of high dispersion is found in the molecular lines 
with increased 0.87 mm continuum emission at a distance of 
2" from the center. The gas here seems to be shocked, however, 
what produces the shock is not clear, an outflow or a spiral arm 
interaction are possible explanations. 

Strong feedback can not be confirmed in the center of 
NGC 1566. Residual velocities that might indicate an outflow 
and which coincide with higher velocity dispersion in molecu¬ 
lar gas. However, these features are aligned with the minor axis 
which can explain the higher dispersion. The LOSV residuals 
might be explained by other means, e.g., streaming motions in¬ 
troduced by the nuclear spiral. These might be signs of di feeding 
of the SMBH. The central ~ 1" harbors possible nuclear star for¬ 
mation as inferred from non-thermal excitation (e.g., UV pho¬ 
tons from young bright stars) on the nucleus. High angular reso¬ 
lution measurements with SINFONI are needed to identify pos¬ 
sible feedback and star formation in the central arcsecond of the 
seemingly waking up Seyfert 1 nucleus. The lack of strong feed¬ 
back could be a reason for the lack of strong star formation in 
the nuclear disk if the nuclear spiral is not able to condense gas 
strong enough to enable star formation. 


In the case of NGC 1566 the distribution and velocities of 
warm and cold molecular gas are very similar. In our first simul¬ 
taneous observation of NUGA sources in the NIR and sub-mm 
in NGC 1433 ( [Smajic et al.[[2Q14] ) the similarities were not as 
high. The distribution was similar but ^^CO(3-2) seemed to be 
a better tracer of the pseudo-ring in NGC 1433. In NGC 1566 
the r = 3" disk is detected in ^^CO(3-2) and H2(1-0)S(1), which 
indicates that the central r ^ 200 pc are warmer than the central 
region of NGC 1433. Additionally, in NGC 1433 the LOSV dif¬ 
fered in value over the full FOV and in distinct regions in LOSV 
and dispersion. This is not the case for NGC 1566 where hot and 
cold gas show similar LOS Vs. 
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Appendix A: Figures 

We have attached the line maps to all important emission lines 
that were detected except for the H 2 Q-branch. The Q-branch is 
situated at the end of K-band and therefore the maps show noise 
features as well. However, the flux and FWHM of these lines 
is given in table for the regions discussed here. Additionally, 


several line ratio maps are presented in Fig. A.l 

Figure [Al2| shows the flux, FWHM, and EW of all other de¬ 
tected H 2 lines detected up to 2.3 //m. The flux distribution is 
very similar in all H 2 lines and peaks on the nucleus. The FWHM 
is increased in the center along the minor axis and H2(1-0)S(3) 
shows an increased width in the northeast and a decreased width 
in the northwest, similar to H2(1-0)S(1). The nuclear spiral is 
detected in all EW maps which show an increased EW at the 
position of the spiral density wave. 


Figure |A.3| shows the LOSV of all H 2 lines from |A.2| and 
additionally the LOSV maps of the detected ionized gas. The 
LOSV in the molecular lines is very similar, as are flux and EW, 
indicating that they all originate from the same region, i.e. the 
nuclear molecular gas disk. For the ionized lines it is hard to 
tell because the region in which they are detected is not large 
enough to give information about the spatial velocity distribu¬ 
tion. However, the LOSV flelds do not look very different from 
the molecular LOSV flelds (e.g., Bry). 
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Fig. A.l - Shown are the logarithmic line ratios of H2(1-0)S(1) over Bry and [Fen] over Bry in (a) and|(b)| Panel (c) shows the logarithmic 
H2(1-0)S(1) over He i ratio. Panel|(d)|(e)|(f)|show H2(2-1)S(1) over H2(1-0)S(1), H2(1-0)S(2) over H2(1-0)S(0) and H2(1-0)S(3) over HzCl- 
O)S(l). 
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Fig. A.2 - Molecular hydrogen emission lines. From left to right: flux [10 W m ^], FWHM (corrected for instrumental broadening) 
[ km s-'l and EW [A] maps of, from top to bottom: H2(1-0)S(3), H2(1-0)S(2), H2(1-0)S(0), and H2(2-1)S(1). Article number, page 19 of 
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Fig. A.3 - Presented are theLOSV [km s '] maps of, from left to right and top to bottom: [Feii], Her, narrow Bry, H2(1-0)S(3), H2(1-0)S(2), 
H2(1-0)S(0), H2(2-1)S(1). 
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